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Using the assumptions of an incompressible and viscous flow at large Reynolds number, 
we derive the evolution equation for the wave action density of an instability wave 
traveling on top of a laminar free-shear flow. The instability is considered to be 
viscous; the purpose of the present work is to include the cumulative effect of the 
(locally) small viscous correction to the wave, over length and time scales on which the 
underlying base flow appears inhomogeneous owing to its viscous diffusion. As such, 
we generalize our previous work for inviscid waves. This generalization appears as an 
additional (but usually non-negligible) term in the equation for the wave action. The 
basic structure of the equation remains unaltered. 


1 . Introduction 

The propagation of modal waves on slightly nonuniform or unsteady base flows can 
produce many interesting effects. Two important entities that describe these effects are 
the wave action -A and its flux G_/£, where G is the group velocity. In the classical 
theories for conservative systems, these two entities are related by a "continuity 
equation" that expresses a conservation law over the long spatio-temporal scales 
associated with the nonuniformity (or inhomogeneity) of the base flow (Whitham 1974, 
Hayes 1970). 

Over the years, there have been many attempts to generalize these concepts to non- 
conservative systems, specifically, to instability waves riding on top of "diverging" base 
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flows. Invariably, these generalizations are formulated in terms of an amplitude 
function A, rather than the wave action jC (see Ho & Huerre 1984 and references cited 
therein). Therefore, even in the simplest case, any correspondence between classical 
kinematic wave theory and the stability of non-parallel base flows is completely lost. 
The equation for the amplitude contains a large number of terms that are difficult to 
interpret, but whose combinations in a certain way form the wave action and the group 
velocity. 

In an attempt to provide such a correspondence, the author formulated the 
linearized instability of non-parallel flows in terms of a complex wave action density, 
jt* and derived an evolution equation for this entity (Balsa 1989; also referred to as (B) 
herein). While the derivation in (B) is very general, it does invoke the assumption that 
the instability wave is inviscid; this assumption turns out to be unnecessary for a large 
class of flows and the analysis of (B) goes through for a viscous instability wave as 
well. Our purpose is to sketch out this new analysis. 

In §2, we describe the underlying base flow (essentially a free-shear layer) that is 
"diverging" or "non-parallel," owing to viscous diffusion. The velocity profile 
associated with this flow is assumed to possess an inflection point; therefore, the 
instability wave is nearly inviscid locally. This will be clear from the equations 
presented in §2; it is the cumulative effect of the small viscous terms that we wish to 
capture, both for the base and perturbed flows. 

In §3, we outline the solutions for the perturbations in terms of a high-frequency 
ansatz. A displacement variable a is used in place of the velocity u. The analysis 
follows that in (B) very closely; we omit most of the details and algebra for the sake of 
brevity. A discussion of these results can be found in §4. 

2. Formulation of the problem 

Consider an incompressible mixing layer for the base flow with velocity U B and 
pressure pg, respectively. The underlying Cartesian coordinate system, x, and the 
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geometry of the problem are illustrated in figure I. All variables are nondimensional: 
x =» (x, y, z) and time, t. are normalized by a characteristic thickness (- L ref ) and a 
transit time (=* L Te{ /U Te{ ) associated with the mixing layer, respectively. Here, £/ ref 
denotes a typical fluid speed (say, that in the upper external stream); this is used to 
nondimensionalize all fluid velocities. As is customary, pressures are measured in 
terms of p£/ ref (p - const. - fluid density). The relevant Reynolds number is Re - 
^refA-ef/" » 1, where v - const, denotes the kinematic viscosity. 

The base flow may be written, to the required order of accuracy, as 

U B - [1 + 0(e 2 )]U(y, {. r) + [1 + 0(£ 2 )]eF(y, |, r)e 2 (la) 

Pg = const. + 0(e 2 ) (16) 

where e = Re* 1 « 1. (ej, e 2 . e 3 ) denote the unit vectors along the coordinate directions, 
and the component of U along the p-axis vanishes. The slow spatio-temporal variables, 
arising from the viscous diffusion of the base flow, are 


£ - €(xej + ze 3 ) , t - a 


(lc) 


It is well known that, under the stated assumptions, the horizontal component of the 
base flow velocity, U, satisfies the boundary layer equations 


where 


du _ a 2 u 

Dt dy 2 


D. 

Dt 




U P 


d_ 

dx 


while the physical law of mass conservation requires 


(2a) 

(2b) 


d_ 

dx 


Up 


= 0 


(2c) 


We shall also use the definition 
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f(y- i. r) = (2d) 

which (i.e., the second derivative of U) plays an important role in the stability analysis. 

Let us assume that this base flow is perturbed by an arbitrary disturbance whose 
(perturbation) velocity and pressure fields are denoted by u - u(x, t) and p - p(x, t), 
respectively. The relevant linearized equations for these quantities are 

Momentum: ^ + u • ^—2- = - + eV 2 u (3a) 

Dr dx dx v ’ 

Continuity: ^ * u - 0 (3d) 

where V 2 denotes the Laplacian operator in (x, y, z) and we have assumed that the 
disturbances are viscous. Our interest is in the case where the perturbations are the 
instability waves of the base flow; quite generally, these waves are nearly inviscid 
when the Reynolds number is larger than about 500 (Betchov & Szewczyk 1963) and the 
base-flow velocity profile possesses an inflection point. 

Since the disturbance equations (3) are nondimensionalized by characteristic length 
and time scales that are representative of instability modes, the viscous term on the 
right-hand side of (3a) is indeed small. This term represents a small (local) correction to 
the instability modes; however, its cumulative effect over the long spatio-temporal 
scales, ( and t, is formally of order unity (Hultgren 1992). We shall incorporate this 
effect into the conservation equation for the wave action density, thereby generalizing 
the result of (B). 


3. Wave action density 

In order to obtain this generalization in the simplest way, we replace the 
disturbance velocity u by a displacement variable a - a(x, t) such that 


u 



9Ub 

dx 


(4) 


A physical interpretation for o - a(x, t) is this: Consider a fluid particle of fixed 
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identity that occupies the point x (at time t) in the unperturbed base flow. The position 
of this same particle in the perturbed flow (at time t) is defined to be (x + a). In other 
words, the particle that would have been at x in the base flow is actually in a slightly 
different position (namely, x + a) in the disturbed flow. 

Using (4) in (3a, b), we obtain 

+ + fV2 u (5a) 

D t 2 dx dx 

£ ■ a - 0 (5b) 

where it turns out to be more convenient to leave the viscous term in its original form. 
The error in (5a) is of 0(e 2 ). It should be noted that (5) is a system of linear partial 
differential equations whose coefficients (i.e., U B in D/Dr and OF) vary slowly with x, 
z, and t. 

Clearly then, in order to obtain an asymptotic solution of (5) when e « 1, we use a 
high-frequency ansatz 

a - a(x, t) - + eaW + ... ]exp (/0/e) + cc (6) 

where / - V~l and cc stands for the complex conjugate of the term immediately 
preceding it. Similar expansions hold for u and p\ 0 - 0({. r) is the complex phase and 
otti) . a(i)(y, {, r) for /' - 0, 1. 

After substituting (6) into (5) and requiring a nontrivial solution for (-)^ with 
vanishing values at y - ±oo. we obtain the evolution equation for the phase 

- - Q(V0, l r) (7a) 

where Q = S2(k, £, t) is, in principle, a known function of its arguments expressing the 
complex local frequency, co - Q ■ -d<j>/d t, as a function of the complex local wave 
vector, k - V0. We call Q, the dispersion relation and note that V denotes the gradient 
operator in propagation space. 
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V " e > | + e 3 | {lb) 

where £ - £ej + fe 3 . Both the frequency and the wave vector, (go, k), are local in the 
sense that they depend on | and t due to the viscous spreading (and possible 
unsteadiness) of the base flow. 

Furthermore, having satisfied the dispersion relation via (7a), the lowest-order 
solutions (ae(°) , //°), u®) must be proportional to an instability mode that can be 
supported by the base flow. If p m = p m (y, r, k) is the pressure mode at wave vector 
k, then we must have 


pM - A(t r) p m (y, $. r, V0) (8a) 

where A - /4(£, t) is the slowly varying complex amplitude of the disturbance. 
Similarly, 


a(°) - /«{. rX0e 2 - iW) . W • e 2 - 0 (8 b) 

where (W, Q) represent the horizontal and vertical components of the displacement 
mode. They are expressible as 

Q * GO'. {.t)> \ (8c) 

<4 dy 

W - W (y.(.r)-- Pm (8rf) 

and p m on the right-hand sides of (8c, d) is to be evaluated at k - V0. The Doppler- 
shifted frequency, co 0 , is given by 

Wq = o>o( y< t) = - + U • V*j (8e) 

The algebra is omitted for brevity; these results are discussed in more detail in (B). At 
this stage of the analysis, the amplitude, A * /4({, r), is an arbitrary function of its 


arguments. 
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In order to obtain the evolution equation for the wave action, we must proceed to 
the next order in e. Here, aO) and //*) will obey linear inhomogeneous equations; the 
forcing terms in these equations arise from the variation of the lowest-order solutions 
on the long spatio-temporal scales, £ and t. In addition, the viscous term in (5a) will 
also contribute to the forcing in the present analysis; this is the only new term that was 
not already discussed in (B). 

Briefly, at 0(e ) , the governing equations for and are 

3p(*) 

- + /V0^*) + e 2 w - terms in (B) + R (9a) 

/V0 • aO) + e 2 • “ terms in (B) (96) 

where R represents the new term that arises from the viscous correction in (5a). It is 
convenient to separate this term into its vertical and horizontal components via 

R * e 2 R± + /Rj| , e 2 * Rj| ■ 0 (9c) 


Employing (4) and (8), we find 


R i - -iA 


dyZ 


V(j> ■ V0 (o> 0 G) 


/R,l - -A 


|5- v * v *](“°w + Qf] 


(10a) 

( 106 ) 


The notation is as follows: The expression "terms in (B)" stands for those forcing 
terms that have been given already by Balsa (1989); we shall not repeat them here 
because they are lengthy and their contribution to the evolution equation for the wave 
action has been calculated in (B). 

Observe that 


^ - V0 • V0 


(ID 


is a representation of the Laplacian operator within the high-frequency ansatz; within 



this same representation, u 0 Q and w 0 W are suitable time derivatives of the fluid 
displacements in cross and propagation spaces. These derivatives provide the lowest- 
order velocity components in these spaces (essentially u^). with the exception that the 
vertical displacement, Q, of a material volume in the nonuniform horizontal base flow, 
U - U(y). also produces a perturbation velocity in propagation space in proportion to 
Q(dU/dy). These effects are clearly contained in (10). 

It is well known that the (*)B) solutions exist if the forcing terms in (9a. b) satisfy a 
solvability condition. When this condition is applied to the "terms in (B)," we obtain 
the result for the wave action density as presented in equation (35) of (B). This 
equation, in our current notation (note that <f> of (B) is written as i<f> herein), becomes 

f + °° 

+ V • (GX) + A 2 \ Q ^ • W dy - 0 (12, also (35) in B) 

97 J_oo dy 

where X and G are defined in (14 b.d). Roughly speaking. (12) contains the effects of 
inviscid instability modes interacting with slightly inhomogeneous base flows. 

On the other hand, when the solvability condition is applied to all the terms on the 
right-hand sides of (9a, b), we obtain the equation for the wave action density as 
presented below (14). This new equation, the main result of this paper, describes the 
effects of slightly viscous instability modes interacting with inhomogeneous base flows. 
Since the solvability condition is linear in the forcing terms, the new equation for the 
wave action contains all the terms presented in (12) plus a new term arising from 
inclusion of R in this condition. This is tantamount to adding the integral (see (31) of 
(B)) 

r 

-iA « QR X - W • R 1l ) dy (13) 

J_oo 

to the left-hand side of (12). After some algebra, the final result may be written as 


%^ + V • (GX) + HX - 0 


(14a) 
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where 


and 


r 

J-oo 


Jt - X(£, t) - A 2 W 0 (G 2 + W • W) dy 


( 146 ) 


r +oo 

V<f> • &{Q 2 + W • W )dy 

J-oo 


H = /ft*, t) 


+ °° 

co 0 «3 2 

J-OO 


+ W • W)dy 


(14c) 


G „ 9Q(k. & r ) 

dk 


k-Vtf 


(14 d) 


We recognize G as the (complex) group velocity. The complex function H 

represents the interaction between the (viscous) wave and the non-uniform base flow; in 
the classical theories, H - 0, thereby reducing (14a) to the familiar conservation 
equation for the wave action, jC (but not for the amplitude A\). 

We emphasize that the integral in (13) contributes only to the interaction function 
H. Also, there are many equivalent ways of writing this contribution (as well as H 
itself; see §4). In order to cast (13) into a simple and useful form, it is necessary to 
recall that (u)q Q) satisfies Rayleigh's equation for the vertical velocity, the pressure 
mode, p m , satisfies Rayleigh's equation for the pressure, and Q and W are expressible 
in terms of p m via (8c, d). The accompanying algebra is reasonably straightforward but, 
perhaps, not entirely trivial. 


4. Discussion and conclusions 

In summary, (14a) provides the evolution equation for the complex wave action 
density, yt, of viscous instability waves interacting with slightly inhomogeneous base 
flows arising from the viscous diffusion of the latter. The complex interaction function, 
H , contains additively both the cumulative effects of the inhomogeneity and those of the 
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viscous correction to the wave. As such, it is the most general result that can be 
obtained within the framework of small-disturbance theory. Note that (14a) is to be 
solved "simultaneously" with the equation for the phase, (7a), in order to provide 30/dr 
and V<f> in the former (though, strictly speaking, these equations are decoupled in the 
linearized theory). 

It is instructive to recast the numerator of (14c) in a form that separates the effects 
of the base flow inhomogeneity from those of the viscous correction for the waves. 
After using some of the tricks mentioned at the end of §3 and performing a couple of 
integrations by parts, we find 


r 

J -oo 


V0 * ^(Q 2 + w • W )dy 

-oo 


-+oo p+OO 

Q • W dy - 2 -j- . Wdy (15) 

J_oo dy J-oo w 0 oy 


where the first term on the right-hand side of (15) is in complete agreement with the 
"inviscid" result in (12). Clearly, the second term is entirely responsible for the viscous 
correction to the instability waves. 

Unfortunately, calculation of the terms in (15) can be done only numerically: it 
entails solving the boundary layer equations (2) for U, solving the inviscid Rayleigh 
equation for the dispersion relation Q - Q(k, £, r) and for the displacement modes Q and 
W, obtaining the phase <)> « $£, r) via (7a), and then evaluating the integrals in (15) 
numerically. This is beyond the scope of the present work. 

Nevertheless, it is possible to analytically calculate the terms in (15), as well as 
those in (14c), for a very special case. In order to do this, we think of H - H(£, t) as a 
function (say, 9T) in the extended space (k, £, r), where 


5T - £, t). 


r 

k • 5»-(Q2 + w • W )dy 

J-oo 

~r 

fi 0 (G 2 + W • W )dy 

J-oo 


( 16 a) 


and the inviscid displacement modes, Q and W, are to be obtained at wave vector k. 
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In addition, 

n 0 - Q - k • U (166) 

so that 

H - H(t r) - 5T(V0. £. r) (16c) 

with the understanding that d<f>/d r - -Q(V0, £. r). Similar remarks hold for (15). 

The special case alluded to previously occurs when, at the point and time of 
interest (which define £ and r). 


U - (tanh y)e l (17 a) 

and the local mode is precisely neutrak. k • k - 1 (Michalke 1964). In this case, we 
have closed-form solutions for the modes and the (singular) integrals given above may 
be evaluated analytically by a well-known contour deformation in the complex y-plane. 
The result is 


9C - 0 (176) 

Evidently, the contributions to 5T arising from the effects of the base flow 
inhomogeneity and those of the viscous correction to the instability mode (as separated 
out in (15)) cancel in this special case; these contributions are +8, respectively. This 
suggests that there may be a (partial) cancellation arising from these two effects, even in 
the general case. Very roughly, we expect the inhomogeneity of the base flow to be 
responsible for an increase in the wave action and the viscous correction to the 
instability to produce a decrease in jC. 

The author is grateful for the support provided by the Institute for Computational 
Mechanics in Propulsion (ICOMP) and the hospitality provided by the NASA Lewis 
Research Center during a summer (1994) visit. Discussions with Dr. L. S. Hultgren 
were extremely helpful in generalizing the author's previous work. 
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Figure I . Geometry of the problem (not to scale), y - cross space; (x, z ) - propagation 
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